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INTRODUCTION 


Photoref ractive  nonlinear  effects  such  as  two-wave  mixing  and 
phase  conjugation  can  be  utilized  for  a  variety  of  image  and 
signal  processing  applications.  The  useful  spectral  range  of 
those  crystals  which  offer  high  gain  is  limited  to  visible 
wavelengths.  Although  photoref ractive  effects  have  been  observed 
in  the  near  infrared  in  BaTiO^  (barium  titanate) ,  the  sensitivity 
has  been  too  low  for  practical  application  at  these  wavelengths. 

A  possible  technique  for  enhancing  the  photoref ractive  effect 
in  BaTiO^  at  laser  diode  wavelengths  was  based  on  the  concept  of 
secondary  photorefractive  centers  introduced  by  Brost  et  al.[l]. 
New  (unpublished)  data  had  indicated  that  the  population  of 
secondary  centers  would  increase  the  absorption  coefficient  at 
800  nm.  The  concept  was  to  populate  these  centers  in  a  manner 
which  could  increase  the  sensitivity. 

This  report  documents  efforts  to  better  understand  the 
characteristics  of  the  secondary  centers  and  to  explore  an 
enhancement  technique  based  on  their  presence  in  BaTiO^.  In 
section  II  a  background  on  photorefraction  and  secondary  centers 
is  given.  Section  III  describes  photo-induced  absorption  and 
two-beam  coupling  experiments  designed  to  yield  information  on 
the  secondary  centers.  In  section  IV  photorefractive  enhancement 
studies,  theoretical  and  experimental,  are  described.  Although 
the  desired  enhancements  were  not  realized,  these  efforts  have 
provided  useful  information  about  the  secondary  centers.  In 
addition  the  theoretical  development  has  applications  in  other 


areas.  In  section  V  additional  photcref ract ive  charactpr ist ics 
of  sublinear  intensity  dependent  rise  time  and  temperature 
dependent  photorefractive  effects  resulting  from  the  presence  of 
secondary  centers  are  discussed. 

II.  BACKGROUND 

The  physical  process  involved  in  the  photorefractive  effect  is 
a  light  induced  redistribution  of  charges  among  traps.  The 
resulting  space  charge  field  then  modulates  the  refractive  index 
through  the  linear  electro-optic  effect.  Based  on  the  simple 
band  transport  model  of  a  single  set  of  recombination  centers 
Kuhktarev  et  al,[2]  derived  solutions  for  formation  of  the 
photorefractive  grating.  These  solutions  have  provided  fairly 
good  agreement  with  many  experiments.  In  an  effort  to  account 
for  intensity  dependent  absorption  and  photorefractive  effects  in 
BaTi03  Brest  et  al.[l]  introduced  the  secondary  center  model. 
Secondary  photorefractive  centers  were  traps  which  were  highly 
thermalized  at  room  temperature  but  became  populated  under  the 
illumination  by  visible  light.  This  model  was  based  on  one  set 
of  secondary  centers.  Brost  et  al.[l]  considered  only  absorption 
effects  at  514  nm.  In  this  effort,  the  secondary  center  model 
for  BaTiO^  was  extended  to  two  species  of  secondary  centers,  as 
suggested  by  the  photo-induced  absorption  data.  This  model  is 
shown  schematically  in  Figure  l.  The  concept  for  enhancing  the 
photorefractive  effect  in  BaTiO^  which  utilizes  the  secondary 
centers  is  shown  in  Figure  2.  The  concept  is  to  illuminate  the 
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Figure  1.  Schematic  diagram  of  secondary  center  model  for 
BaTiO^ . 
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crystal  with  pump  beam  of  visible  wavelength  while  writing  the 
photoref ract ive  grating  with  a  cw  laser  at  a  near  infrared 
wavelength.  It  was  hoped  that  the  population  of  the  secondary 
centers  by  the  pump  beam  could  result  in  an  enhancement  through 
an  increased  sensitivity  associated  with  the  increased  absorption 
in  the  near  infrared. 

The  response  of  the  photoref ractive  material  to  a  spatially 
modulated  light  intensity  pattern  is  well  described  by  a  set  of 
nonlinear  differential  equations.  The  relevant  equations  for  the 
secondary  center  shown  in  Figure  1  are 


dNj^/dt  -  -  N^)N)^  (1) 

dN2/dt  =  ”(52111  +S22I2  ■*'  ^2^^2  >2^^2T  ~  ^2^^h 

dN3/dt  =  “(Siili  +  522^2  ^3)^3  +  I3(^3T  ~  ^3^^h 

dNj^/dt  =  ~dNi/dt  -"dN2/dt  -dN3/dt  (4) 

aNj^/dt  =  dNj^/dt  -  V-jj^  (5} 

ih  =■  <^> 

V*E  =  p/€fQ,  (7) 


where  Nix  is  the  total  number  density  of  photore ^ract i ve 
centers,  N21.  and  N3rp  are  the  total  number  densities  of  the 
secondary  centers,  is  the  number  density  of  unionized  charge 
donors,  is  the  number  density  of  free  carriers  (these 

equations  are  written  for  positively  charged  free  carriers) , 
is  the  free  carrier  current  density,  Ii  and  I2  are  the  light 
intensities  of  two  different  wavelengths,  2  the  photon- 
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absorption  cross  sections  for  light  of  wavelengths  1  and  2,  is 
the  recombination  coefficient,  is  the  thermal  ionization  rate, 
is  the  carrier  mobility,  D  is  the  diffusion  coefficient  given 
by  D  =  ^^kgT/q,  E  is  the  electric  field,  t  is  the  dielectric 
constant,  is  the  permittivity  of  free  space,  p  is  the  charge 
density.  These  equations  are  subject  to  the  constraint  of 
charge  conservation 

Ni+  N2+  N3  +  Nj^  -  Nq  =  0,  (8) 

where  Nq  =  +  N2fj  '^^3d'^  ^hd  ^hd'  ^id 

initial  dark  concentration  levels  before  the  light  was  turned  on. 
The  two  beam  coupling  gain  coefficient  is  then  related  to  the 
photorefractive  space  charge  field  through 

r  =  n^mr^fflmCEj)  (9) 

where  n  is  the  refractive  index,  r^^^  is  the  effective  electro¬ 
optic  coefficient  and  is  the  fundamental  component  of  the 
space  charge  field. 

For  photorefractive  effects,  the  light  intensity  pattern  is 
assumed  to  be  spatially  periodic  and  of  the  following  form: 

I(x,t)  =  +  msin[Kx]},  (10) 


where  is  the  total  intensity,  m  is  the  modulation  index,  K  is 
the  grating  wave  number.  For  two  vertically  polarized  input 
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waves  of  intensities  Ig  and  Ipj,  m  -  (B  +  l)  ,  where  5  =  Ip^/Ig. 

Examination  of  equations  (I-IO)  indicate  t]\at  there  are 
eighteen  input  parameters,  eight  of  which  are  directly  related  to 
the  secondary  centers.  Most  of  these  parameters  associated  with 
the  secondary  centers  were  unknown.  Consequently  experiments  were 
conducted  to  improve  the  understanding  of  the  secondary  centers. 

III.  SECONDARY  CENTER  CHARACTERIZATION 

This  section  describes  results  of  experiments  designed  to 
further  characterize  the  secondary  centers  in  BaTiO^.  By 
analyzing  many  separate,  but  related  experiments  and  improved 
picture  of  the  secondary  centers  could  be  obtained.  These 
included  photo- i nduced  absorption  and  two-beam  coupling 
experiments . 

Photo-induced  absorption 

A  important  characteristic  of  the  secondary  centers  is  the 
presence  of  photo-induced  absorption.  Motes  and  Kim  [3]  first 
reported  the  existence  of  what  was  called  intensity  dependent 
absorption.  In  their  experiments,  the  absorption  coefficient  of 
an  argon  ion  laser  was  found  to  depend  on  its  intensity.  Brost 
et  al.[l]  later  explained  this  absorption  in  terms  of  the 
secondary  center  model.  In  this  report,  the  intensity  dependent 
experiments  were  extended  to  measure  the  photo-induced  absorption 
spectrum. 

A  schematic  diagram  of  the  experiment  is  shown  in  Figure  3. 
Photo-induced  absorption  was  generated  in  the  barium  titanate 
crystal  by  illumination  with  the  514  nm  line  of  an  argon-ion 
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laser.  The  transmitted  power  of  a  weak  signal  beam  (I3)  was 
monitored  by  a  large  area  photodiode.  The  signal  beam  was 
obtained  from  a  monochromator  illuminator  consisting  of  a 
tungsten  halogen  lamp  and  a  grating  monochromator.  It  was 
focused  to  approximately  1.5  mm  diameter  in  the  barium  titanate 
crystal.  The  intensity  of  the  signal  beam  varied  between  3  and 
100  W/cm^,  depending  upon  wavelength  and  polarization.  The  pump 
beam  was  expanded  to  fully  illuminate  the  face  of  the  crystal, 
which  was  5x5  mm^ .  Care  was  taken  to  minimize  scattered  light 
from  the  pump  beam  reaching  the  photodiode  by  use  of  filters  and 
beam  stops. 

The  magnitude  of  Ig  was  reduced  in  the  presence  of  the  argon 
laser  pump  beam  (Ip)  due  to  the  photo-induced  absorption.  The 
photo-induced  component  of  the  absorption  coefficient,  Uj,  was 
determined  using  the  relation 

Ir.  (with  Ip  on)  =  exp(-  a^-d)  ,  (11) 

I3  (with  Ip  off) 

where  d  is  the  interaction  length.  A  typical  spectrum  of  the 
induced  steady-state  absorption  is  shown  in  Figure  4.  The 
induced  absorption  spectrum  did  not  depend  on  the  polarization 
of  the  pump  light.  An  anisotropy  with  respect  to  the 
polarization  of  the  signal  light  was  observed.  The  magnitude  of 
the  induced  absorption  for  extraordinary  polarized  light  was 
about  35%  less  than  chat  for  ordinary  polarization.  A  germanium 
detector  was  used  to  check  the  spectrum  between  l.O  and  1.8  m. 
No  photo-induced  absorption  peaks  were  found  in  that  region. 

The  temporal  (rise  and  decay)  characteristics  of  the  induced 
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absorption  were  measured  by  turning  the  pump  beam  on  for  a  time 
period  sufficient  to  establish  steady-state  conditions  while 
monitoring  the  signal  beam  transmission,  as  shown  in  Figures  5 
and  6.  The  rise  time  for  the  absorption  to  reach  steady-state 
depended  upon  the  pump  intensity,  and  was  typically  about  an 
order  of  magnitude  faster  than  the  photoref ract i ve  grating 
formation  time.  The  dark  decay  characteristics  depended  upon 
the  signal  beam  wavelength.  At  440  nm  the  induced  absorption  was 
characterized  by  a  fast  decay  rate  of  approximately  60  ms  time 
constant.  At  800  nm  two  distinct  decay  rates  were  apparent:  a 
fast  decay  similar  to  that  at  440  nm,  and  a  slow  component,  with 
a  time  constant  of  about  6  s.  The  spectrum  could  thus  be 
separated  into  two  components,  "fast”  and  "slow",  as  shown  in 
Figure  7.  Here  the  slow  decay  spectrum  was  defined  as  the 
absorption  remaining  2  seconds  after  the  pump  beam  was  turned 
off,  and  the  fast  decay  was  the  remainder  of  the  steady  state 
absorption.  The  broad  infrared  steady-state  absorption  peak 
could  then  be  separated  into  two  absorption  peaks  at  775  and  820 
nm. 

The  relative  magnitudes  of  the  "slow"  and  "fast"  decay 
components  of  the  induced  absorption  depended  upon  the  intensity 
of  the  pump  beam.  This  is  seen  in  the  data  shown  in  Figure  8 
for  a  800  nm  wavelength  signal  beam  and  a  pump  intensity  of  4.3 
mW/cm^.  Only  the  "slow"  absorption  component  was  induced  with 
this  low  pump  intensity.  The  magnitude  of  the  slow  component 
saturated  at  a  pump  intensity  of  about  10  mW/cm^.  The  fast  decay 
component  required  higher  laser  intensities  and  was  observed  only 
for  pump  intensities  above  30  mW/cm^,  and  saturated  at  about  10 
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Figure  6.  Signal  beam  transmission  (Ig)  versus  time  for  800  nm 
wavelength  signal. 
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Figure  7.  Time-resolved  photo-induced  absorption  spectrum  in 
BaTiOo . 
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Figure  8.  Signal  beam  transmission  (Ig)  versus  time  for  800  nm 
wavelength  signal  with  low  intensity  pump  beam. 


W/cm^.  The  intensity  dependence  is  shown  in  Figure  9. 

This  photo-induced  absorption  data  indicates  that  there  are 
two  dominant  secondary  photorefractive  centers.  The  magnitude  of 
the  photo  induced  absorption  coefficient  at  saturation  provides 
useful  information  about  the  secondary  centers  through  the 
relation 

“i  =  ^i^iT-  (^2) 

Temperature  dependent  measurements  of  the  photo- induced 
absorption  were  also  preformed.  The  goal  of  these  experiments 
was  to  measure  the  thermal  activation  energy  of  the  secondary 
centers  by  measuring  the  temperature  dependence  of  the  decay  rate 
of  the  photo-induced  absorption.  For  these  experiments  the 
BaTi03  crystal  was  mounted  inside  a  glass  cell  containing 
silicon  oil.  The  cell  was  heated  with  a  hot  plate  and  the 
temperature  was  monitored  with  a  thermo-couple  probe.  Based  on 
the  photo-induced  absorption  spectrum  measurements,  it  was 
possible  to  selectively  probe  the  fast  and  slow  decay  components 
of  the  photo-induced  absorption  corresponding  to  N2  and  N3,  by 
the  choice  of  signal  wavelength.  A  457  nm  line  from  an  argon  ion 
laser  was  used  as  a  probe  of  the  fast  decay  component.  For  the 
slow  decay  component  a  750  nm  diode  laser  was  used. 

The  magnitude  of  the  photo-induced  absorption  decreased  with 
increasing  temperature.  This  is  demonstrated  in  Figure  10  where 
the  temperature  dependence  of  the  photo-induced  absorption 
coefficient  is  plotted  for  a  pump  intensity  of  10  W/cm  .  This 
decrease  in  absorption  is  attributed  to  the  increased  thermal 
ionization  rate  with  increasing  temperature. 
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Figure  9.  Intensity  dependence  of  photo-induces  absorption  at 
various  signal  wavelengths  for  a  514  nm  pump  beam. 
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Figure  10.  Temperature  dependece  of  the  phoro-induced  absorp 
coefficient  at  457  nm  for  a  10  W/cm^  pump  intensity  at  514  nm 


Figure  11.  Photo- induced  absorption  decay  time  as  functi 
temperature. 


The  results  of  the  thermal  decay  experiments  are  shown  in 
Figure  11.  Here  the  thermal  decay  time  is  plotted  on  an  Ahrenius 
plot  as  the  natural  log  of  1/ft  versus  1/T.  The  linear 
relationship  on  this  plot  indicates  that  a  thermal  activation 
energy  in  the  form  of  equation  11.  From  this  data  the  thermal 
activation  energies  for  the  fast  and  slow  decay  components  were 
determined  to  be  0.5  and  0.8  eV,  respectively. 

Two  Beam  Coupling  Experiments 

Additional  information  about  the  secondary  centers  can  be 
obtained  from  two  beam  coupling  experiments.  An  estimate  of  the 
total  concentrations  N2ip  and  N^-j,  can  be  obtained  from  the 
intensity-dependent  two  beam  coupling  data.  According  to  the 
band  transport  model  of  photorefraction  the  maximum  two-wave 
mixing  gain  coefficient  will  occur  at  a  grating  period  that  is 
equal  to  the  Debye  screening  length  Lg,  which  is  given  by 

Lg2  =  ee^kgT/q^Ng.  (13) 

The  effective  empty  trap  concentration  Ng  can  then  be  determined 
from  a  photorefractive  measurement  of  Lg.  As  shown  by  Brest  et 
al.  [1],  the  intensity  dependent  population  of  the  secondary 
centers  results  in  an  intensity-dependent  effective  empty  trap 
concentration.  The  difference  in  the  high  and  low  intensity 
values  of  Ng  is  approximately  equal  to  the  secondary  center 
concentration.  The  intensity-dependent  two-beam  coupling  data 
for  a  wavelength  of  514  nm  reported  by  Brest  et  al.  [1]  is  shown 
in  Figure  12.  This  data  indicated  that  the  total  secondary 
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Figure  12.  Intensity  dependent  two  beam  coupling  gain  versus 
grating  period  for  514  nm  wavelength. 
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Figure  13.  Intensity  dependent  two  beam  coupling  gain  versus 
grating  period  for  800  nm  wavelength. 
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center  concentration  was  approximately  1.5  x  cm 

coupling  data  taken  at  800  nm  wavelength  is  shown  in  Figure  13. 
These  results  indicate  that  at  infrared  wavelengths  at  very  low 
intensities  (l  mW/cm  )  secondary  centers  do  not  contribute  to 
the  photoref ract ive  effect,  and  is  approximately  9x10^^  cm~^. 

But  at  somewhat  higher  intensities  (100  mW/cm^  )  the  secondary 
centers  (  corresponding  to  do  contribute  to  the 

photoref ractive  effect.  From  this  data  we  can  estimate  that  the 
secondary  center  concentration  is  about  2.5x10^^  cm~^.  It  is 
also  apparent  that  at  these  moderate  intensity  levels  external 
pumping  of  the  secondary  centers  will  not  result  in  strong 
enhancements  as  the  secondary  centers  are  already  contributing  to 
the  photoref ractive  effect. 

IV.  PHOTOREPRACTIVE  ENHANCEMENTS 

The  photoref ract ive  enhancements  from  utilization  of  the 
secondary  centers  were  explored  theoretically  and  experimentally. 
Although  the  desired  enhancements  were  not  realized,  these 
efforts  were  are  described  below. 

Theory 

The  equations  (1-10)  were  studied  theoretically  to  try  to 
determine  the  optimum  conditions  for  photoref ractive  enhancement. 
Two  approaches  were  developed  to  solve  nonlinearly  coupled 
equations  describing  the  photoref ract ive  effect.  The  first 
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approach  based  on  the  small  m  approximation  was  appropriate  for 
conditions  of  cw  illumination.  To  study  the  effects  of  pulsed 
illumination,  a  second  technique  was  developed  based  on  a  finite 
difference  technique. 

In  the  first  approach  the  small  modulation  approximation  was 
made  in  which  only  the  first  order  terms  are  kept  in  the 
solution.  Solution  of  the  material  equations  proceeded  by 
setting  N^,  N2,  N3,  and  E  to  the  form  X  =  Xq  +  and  solving 
for  the  lowest  Fourier  component  in  the  quasi-steady 
approximation.  The  X^  quantities  are  assumed  to  be  independent  of 
space  and  time,  and  used  as  input  parameters  for  the  equations 
involving  X^^.  This  approach  linearized  the  equations.  The  rate 
equations  were  solved  numerically  to  determine  the  zeroth  order 
terms.  After  eliminating  the  current,  a  set  of  linear  equations 
in  five  variables  remained.  Analytical  expressions  were 
intractable,  and  consequently  the  solutions  were  obtained 
numerical."' y.  This  approach  was  appropriate  only  for  continuous 
illumination.  Photoref ract ive  space  charge  fields  were 
calculated  for  the  conditions  in  which  the  laser  write  beams 
(corresponding  to  Ij^)  were  assumed  to  be  at  800  nm,  while  the 
crystal  was  illuminated  with  a  pump  argon  beam  at  514  nm. 
Parameters  for  the  secondary  centers  were  based  on  measurements 
reported  elsewhere  in  this  report.  A  range  of  laser  intensities 
of  both  the  write  beams  and  the  pump  beams  were  tried.  In  all 
cases,  no  conditions  of  enhancement  could  be  found.  The  grating 
erasing  effect  of  the  argon  pump  beam  was  always  greater  than  any 
enhancement  effects  due  to  the  population  of  the  secondary 
centers. 
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In  order  to  determine  the  characteristics  of  the  space-charge 
field  under  conditions  of  pulsed  illumination  of  the  pump  beam  a 
finite  difference  scheme  was  developed  to  model  the 
photorefractive  grating  formations.  Unfortunately,  this  approach 
was  very  computer  intensive.  The  nonlinearly  coupled  material 
equations  were  very  unstable,  requiring  a  small  time  step.  It 
turned  out  to  be  impractical  to  use  this  approach  model 
interactions  in  BaTi03,  due  to  the  extensive  cpu  time  required. 
However,  this  approach  has  been  successfully  used  to  model 
photorefractive  grating  formations  in  other  crystals. 

In  the  second  approach,  the  microscopic  space-charge  field 
was  calculated  under  the  assumptions  that  the  spatial  symmetry  of 
the  grating  formation  was  one  dimensional  and  periodic  in  the 
grating  vector  direction.  With  these  assumptions,  it  was  only 
necessary  to  calculate  the  charge  redistribution  over  one  grating 
period  of  length,  which  was  divided  into  N  cells.  The  electric 
field,  and  free  carrier,  carrier  current,  and  charge  densities 
were  determined  for  each  cell.  The  current  was  allowed  to  flow 
for  a  time  period  ^t  and  the  change  in  the  free  carrier 
concentration  was  determined  by  the  continuity  equation.  The 
distribution  of  charge  among  free  carriers  and  traps  was  then 
recalculated  and  the  other  quantities  updated. 

The  distribution  of  charges  among  the  traps  and  free  carrier 
concentrations  in  each  cell  were  determined  from  the  rate  Eqs 
(1-3)  carriers  in  the  ith  cell  due  to  the  nonuniform  current 
flow.  Since  changes  at  each  time  step,  the  constraint  given 

by  Eg.  (8)  also  varies  at  each  time  iteration.  Two  techniques 
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were  used  to  solve  the  rate  equations.  The  first  technique  was 
based  on  the  quasi  static  approximation  in  which  the  time  for 
establishing  free  carrier  equilibrium  was  assumed  to  be  much 
faster  than  the  grating  formation  time.  In  this  approximation 
the  rate  equations  were  solved  for  the  steady  state  solution  at 
each  time  step.  This  approach  was  not  cpu  intensive  and  worked 
well  for  cw  illumination,  but  could  not  work  with  pulsed 
illumination.  For  pulsed  illumination,  a  fourth  order  Runge 
Kutta  technique  was  used  to  integrate  the  rate  equations. 

The  free  carrier  current  was  then  determined  for  each  cell 
from  Eq.  (4).  The  photoinduced  space  charge  field  could  be 
determined  by  combining  the  one  dimensional  continuity  equation 
with  Poisson's  equation  to  give 

‘t 

Esc(x,t)  =  -q/«o  G(t),  (14) 

Jo 

where  G(t)  is  determined  from  the  boundary  condition  of  the 
constraint  of  a  constant  applied  voltage  given  by 

’L 

-  E(X,t)dx  =  V,  (15) 

,  0 

Where  E(x,t)  =  Eg^^(x,t)  +  (x,t)  -  V(t)/L  is  the  total 

electrostatic  field  and  L  is  the  crystal  length.  With  the 
assumption  of  periodicity  this  means  that 

'x*+A 

Eg(,(x,t)dx  =  0.  (16) 

x* 

•/ 
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The  space  charge  field  was  determined  by  applying  Eqs.  (14-16) . 

The  one  dimensional  continuity  equation  was  used  in  the 
following  form  to  determine  the  change  in  the  free  carrier 
concentration  in  the  ith  cell: 

=  -(V.j^.)<yt.  (17) 

The  results  from  Eq .  (17)  were  used  to  update  the  charge 
conservation  constraint  given  by  Eq.  (8)  in  the  solution  of  the 
rate  equations  on  the  next  iteration. 

These  calculations  provided  a  temporal  evolution  of  the  space 
charge  field  as  well  as  other  quantities  such  as  free  carrier  and 
charge  densities.  For  stability,  the  time  step  at  each  iteration 
was  required  to  be  less  than  the  free  carrier  lifetime  and 
characteristic  transport  times.  For  parameters  appropriate  to 
BaTi03,  this  meant  time  steps  of  less  than  l  nsec  even  for  cw 
illumination.  When  considering  that  grating  formation  times  can 
be  many  seconds,  this  approach  of  integrating  the  rate  equations 
required  too  much  cpu  time  to  be  of  practical  use.  Attempts  were 
made  to  combine  the  quasi-static  approach  with  the  Runge  Kutta 
integration  to  study  the  effects  of  pulsed  illumination  from  the 
pump  beam  while  writing  the  grating  with  cw  beams,  and  thereby 
reduce  the  cpu  time.  Unfortunately  these  attempts  were 
unsuccessful. 
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Experiment 


Experiments  were  conducted  to  test  out  the  enhancement  concept 
depicted  in  Figure  2.  In  these  experiments  a  800  nm  line  from  a 
Titsaphire  laser  was  used  for  the  two  beam  coupling.  The  signal 
beam  I],  was  monitored  for  enhancements  in  either  the  speed  of  the 
grating  formation  or  the  magnitude  of  the  two-beam  coupling  gain. 
Either  a  doubled  Nd:YAG  laser  or  the  514  nm  line  from  an  argon- 
ion  laser  was  used  as  the  pump  beam. 

In  the  first  set  of  experiments  a  cw  pump  beam  was  used. 
Intensity  levels  on  the  order  of  a  few  mW/cm^  were  used.  As 
expected  from  the  theoretical  modeling  no  enhancements  were 
observed.  The  presence  of  a  cw  pump  beam  had  the  effect  of 
simply  increasing  the  dark  current  and  thereby  lowering  the  gain. 
The  erasure  due  to  the  pump  beam  dominated  any  enhancements 
effects  that  might  be  present. 

The  next  approach  to  enhancement  was  to  pulse  the  pump  beam. 
It  was  anticipated  that  if  the  secondary  centers  could  be 
populated  with  a  short  pulse  of  green  light  then  the  erasure 
effects  would  be  minimized.  It  was  known  from  the  photo-induced 
absorption  experiments  that  the  time  required  to  populate  the 
secondary  centers  was  much  shorter  than  the  time  required  to 
write  the  grating.  A  number  of  pulse  illumination  configurations 
were  tested.  A  mechanical  shutter  was  tested  for  illuminations  > 
greater  than  1  msec.  An  acousto-optic  bragg  cell  was  used  as  a 
shutter  for  illuminations  down  to  a  few  microseconds.  These 
pulsed  illumination  experiments  were  tested  for  a  variety  of 
repetition  rates  and  illumination  intensities.  Although  the 
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erasure  effects  were  much  less  noticeable  than  with  the  cw 
illumination,  no  enhancemets  were  observed. 

It  was  thought  that  if  the  pulsed  illuminations  could  be  made 
sufficiently  short  then  there  would  be  little  or  no  erasure  from 
the  free  carriers.  The  free  carrier  lifetime  in  BaTi03  is 
thought  to  be  on  the  order  of  1  nsec.  The  diffusion  lengths  of 
a  free  carrier  in  BaTiOj  is  approximately  0.02  nm,  which  is  much 
less  than  the  typical  grating  period  of  0.5  urn.  Consequently, 
for  very  short  pulse  illuminations  the  free  carriers  would  be 
trapped  before  traveling  the  distance  of  1/2  grating  period 
needed  to  erase  the  grating.  To  test  the  enhancement  with  short 
pulses  a  Q-switched  Nd:YAG  laser  was  used.  The  pulse  widths  were 
approximately  5  nsec.  Pump  repetition  rates  ranged  from  0.5  to 
10  Hz.  Fluence  levels  of  up  to  50  mJ/cm^  were  used.  Again,  no 
enhancements  in  the  photorefractive  effects  were  observed. 

V.  OTHER  SECONDARY  CENTER  EFFECTS 

In  the  course  of  this  effort  it  became  apparent  that  the 
secondary  photorefractive  centers  were  responsible  for  many 
photorefractive  characteristics  of  BaTi03  which  are  not  explained 
by  the  usual  model  for  photorefraction.  Two  of  these  effects 
are  the  sublinear  intensity  dependence  of  the  grating  formation 
time  constant  and  the  temperature-dependent  two  beam  coupling 
effective  trap  concentration. 
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Sublinear  intensity  dependence 

An  important  issue  in  the  application  of  the  photorefractive 
effect  is  the  speed  of  the  grating  formation.  According  to  the 
widely  used  results  of  Kukhtarev [ 2 ] ,  the  grating  formation  time 
is  expected  to  be  inversely  proportional  to  the  laser  intensity. 
However,  experiments  performed  on  barium  titanate  have  shown  that 
the  response  time  varies  with  intensity  as  where  the 

exponent  x  ranges  from  0.5  to  1.0. [4-6]  In  this  section  we 
show  that  this  sublinear  dependence  of  the  photorefractive 
grating  formation  time  is  explained  by  the  presence  of  secondary 
photorefractive  centers.  For  convenience,  the  derivation  is 
limited  to  the  presence  of  only  one  type  of  secondary  center. 

According  to  the  Kukhtarev  solutions,  the  photorefractive 
response  time  r  can  be  written  as [7]. 


T 


~  ^die 


1  + 

1  +  K^Lg^ 


(18) 


where  is  the  dielectric  relaxation  time  given  by 

^die  =  «€o/‘3^hNhO'  (19) 

Lg  and  Ljj  are  the  screening  and  diffusion  lengths,  respectively, 
given  by 

=  eCokgT/q^Nj,  (2  0) 

and 

=  ^hW"?'  (21) 

K  is  the  grating  wave  number,  Ng  is  the  effective  empty  trap 
density,  is  the  hole  mobility,  Nj^q  is  the  mean  free  carrier 
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concentration,  Tj^  is  the  free  carrier  lifetime,  e  is  the  dielec¬ 
tric  constant.  In  the  single  donor-trap  model,  is  the  only 

intensity-dependent  parameter  which  affects  equation  (18).  It 
is  simply  proportional  to  the  laser  intensity,  and  consequently, 
the  response  time  is  expected  to  vary  inversely  with  intensity - 
The  important  feature  of  the  secondary  center  model  is  the 
finite  thermal  ionization  rate  of  the  secondary  centers  (^2  ~  30 
sec”^)  which  influences  the  free  carrier  and  occupied  trap 
concentrations.  The  photorefractive  response  time  predicted  by 
the  secondary  center  model  can  be  determined  in  a  manner  similar 
to  that  used  by  Valley[8].  We  start  with  the  equations  used 
to  calculate  the  steady-state  space-charge  field: 


dN^/dt  =  -  Nj^)Nh  (22) 

dNj/dt  =  -(S21I1  +S22I2  +  3'2^^2T  “  ^2)^h  (^3) 

dNjj/dt  =  -dNj^/dt  -dNj/dt  (24) 

aN^/dt  =  dNf^/dt  -  V.jj^  (25) 

3h  =  ^h^h®  -  (26) 

V*B  =  p/ecQ,  (27) 


Solution  of  equations  (22  -27)  proceeds  by  setting  Nj^,  N2,  Nj^, 
and  E  to  the  form  X  =  Xg  +  and  solving  for  the  lowest  Fourier 
component  in  the  quasi-steady  approximation[ 1] .  The  Xg  quantities 
are  assumed  to  be  independent  of  space  and  time  while  the  first 
order  quantities  are  assumed  to  be  of  the  form  exp(iKz-rt) .  This 
yields  the  following  set  of  simultaneous  equations  in  the  first 
order  variables: 
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0 


(ri~r)N3^3^  -  ^R1%1  - 
(r2“r)N2i  -  rj^2*^hi  ”  "®2^o”'^20 

/^h^hoAi  -rNii  -rN2i  +  (rD-r)Nhi  = 
(“€eQ/q)Aj^  +  +  N2]^  +  =  0 


(28) 

(29) 

(30) 

(31) 


where  -  S2lo+fii+  i^ho'  ^2  “  ^2^0'^®2'‘'  2^h0'  ^R1  ~ 

^R2  =  2(N2T-N2o)'  =  K^^^kgT/q,  and 

intensity  modulation  index,  Aj^  =  iKE^^,  and  Eq=0.  Nj^q,  N2o»  and 
Nj^q  are  constants  determined  from  the  zeroth  order  equations. 
The  photorefractive  response  time  1/r  is  found  by  setting  the 
determinant  of  left  hand  side  of  equations  (28-31)  equal  to  zero. 
With  the  approximation  I’d  »  r^j_g-r,  this  yields  equation  for  the 
photorefractive  response  time  1/r 

(rD+rR2+rRi)r2  -  [  rg(r3^+r2)  +  rjjie(rj^j^+rjj2)  '^^i^ri  ■‘■^2^ri 

4.  [rgrir2  +  rdierirg2  +  =  o.  (32) 

The  intensity  dependence  of  the  response  time  is  shown  in 
Figure  14.  The  lines  are  theoretical  results  calculated  from 
equation  (32)  and  the  circles  and  triangles  are  experimental 
data.  This  data  was  taken  on  an  as-grown  p-type  barium  titanate 
crystal  using  the  514.5  nm  line  of  an  argon-ion  laser.  The 
grating  wavevector  was  parallel  to  the  C-axis  and  the  light  had 
ordinary  polarization.  The  ratio  of  the  signal  beam  to  the  pump 
beam  was  0.03  to  1.  We  monitored  the  change  in  power  of  the 
signal  beam  in  the  presence  of  the  pump  beam  and  found  the  time 
dependence  of  the  grating  formation  to  be  described  as  a  single 
exponential.  The  power  law  dependencies  determined  from  a  fit 
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INTENSITY  (¥/cm2) 


Figure  14.  Photorefractive  rise  time  versus  incident  laser 
intensity  for  0.4  and  1.2  grating  periods.  The  circles  and 
triangles  are  experimental  data  and  the  lines  are  secondary- 
center  model  theory. 
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the  data  to  an  form  were  for  a  1.2  fj.m  grating  period 

and  i~0.56  jQj.  ^  0.40  nm  grating  period.  The  parameter  values 
used  in  the  calculations  are  listed  below  and  are  the  same  that 
were  used  in  reference  l  to  model  the  steady-state 
characteristics  of  the  same  barium  titanate  crystal.  The  hole 
mobility  was  the  only  new  parameter.  Changing  the  magnitude  of 
the  mobility  did  not  change  the  shape  of  the  curve,  but  did 
determine  the  magnitude  of  the  response  time.  We  found  that  a 
value  of  about  0.25  cm^/V-sec  best  matched  the  magnitude  of  the 
experimental  response  times. 

The  theoretical  curves  were  only  approximately  of  the  form  l“^. 
The  response  times  varied  as  l/I  at  the  low  and  high  intensity 
regimes  but  were  sublinear  at  intermediate  intensities,  the 
region  of  most  experiments.  A  fit  of  the  theoretical  curves 
to  an  functional  form  over  an  intensity  range  of  0.01  to  10.0 
W/cm^  gave  power  law  dependencies  of  i"0*'78  and 
1.2  and  0.4  nm  grating  periods,  respectively. 

It  is  convenient  to  discuss  the  effect  of  the  secondary 
centers  on  the  response  time  by  using  equation  (18).  This 
expression  remains  meaningful  for  the  secondary  center  model;  it 
is  the  values  of  the  dielectric  relaxation  time  and  diffusion 
and  screening  length  which  are  model  dependent.  The  same  results 
were  obtained  using  either  equations  (18)  or  (32) . 

The  free  carrier  concentration  influences  the  photoref ractive 
response  time  through  the  1/Nj^O  dependence  of  the  dielectric 
relaxation  time.  One  important  consequence  of  the  secondary 
center  model  is  the  sublinear  dependence  of  the  free  carrier 
concentration  on  laser  intensity.  This  is  shown  in  Fig.  15.  Here 
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Figure  15.  Free  carrier  concentration  versus  laser  intensity  for 
the  secondary-center  and  single-impurity  models. 


we  plot  the  calculated  steady-state  free  carrier  concentration 
under  uniform  illumination  as  a  function  of  laser  intensity  for 
both  the  single-impurity  and  the  secondary-center  models. 
While  the  free  carrier  concentration  increases  linearly  with 
intensity  for  the  single-impurity  model,  the  presence  of 
secondary  centers  results  in  a  sublinear  dependence. 

The  response  time  also  depends  on  the  screening  and  diffusion 
lengths,  as  seen  by  examination  of  equation  (18).  From  this 
equation  we  see  that  at  long  grating  periods  (>  4-0  in  BaTiO^) 
the  response  time  is  approximately  equal  to  the  dielectric 
relaxation  time.  At  smaller  grating  periods  the  second  term  is 
significant.  In  BaTiO^  the  screening  length  is  typically  larger 
than  the  diffusion  length  and,  therefore,  the  response  time 
decreases  as  the  grating  period  decreases. 

The  diffusion  length  may  vary  as  a  result  of  an  intensity 
dependence  in  the  free  carrier  lifetime,  due  to  changes  in  the 
occupied  trap  concentrations.  This  is  shown  in  Fig.  16  where  we 
plot  the  lifetime  versus  intensity.  A  more  important  effect  is 
the  intensity  dependence  of  the  screening  length.  According  to 
equation  (20) ,  Lg^  is  inversely  proportional  to  Ng.  As  was  shown 
in  reference  1,  the  effective  empty  trap  concentration  increased 
with  laser  intensity  due  to  the  contribution  of  the  secondary 
centers  to  the  redistribution  of  charges.  This  effect  is  shown 
in  Fig.  16  where  we  plot  Ng  versus  laser  intensity.  The 
effective  empty  trap  concentration  increased  by  about  a  factor  of 
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Figure  16.  Effective  empty  trap  concentration  and  free  carrier 
lifetime  versus  laser  incensity  for  the  secondary  center  model. 


In  Figure  17  we  plot  the  calculated  power  law  dependence  of 
the  photorefractive  rise  time  versus  intensity  as  a  function  of 
grating  period.  This  was  determined  by  fitting  the  theoretical 
curves  to  an  1“^  functional  form  for  intensities  between  0.01  and 
10  W/cm^.  At  long  grating  periods  the  power  law  dependence  is 
the  same  as  that  of  the  dielectric  relaxation  time  and  is 
determined  by  the  sublinear  free  carrier  dependence,  as  shown  in 
Fig  15.  At  small  grating  periods  the  exponent  decreases  with 
decreasing  grating  period.  This  is  a  result  of  the  intensity- 
dependence  of  the  characteristic  lengths,  primarily  that  of  the 
screening  length  due  to  the  intensity-dependent  effective  empty 
trap  concentration. 

Parameter  values  used  in  the  above  calculations  are:  = 
2.5X10“^^  cm^,  S2  =  5X10"^®  cm^,  ^  =  5X10"®  cm®s”^,  ^  = 
2.5X10“®  cm®sec“^,  6^^  =  1.5X10“^  sec"^,  ^2  =  sec"^, 
N;ILT  =  2.6X10^®cm“®,  N3^(0)  =  2.57X10^®  cm“®,  N2T  =  1.5X10^’^  cm  % 
N2(0)  =  3.8X10^®  cm“®,  Nj^(O)  =  3X10^  cm“® ,  e  =  168,  and  = 
0.25  cm^V'^s”^. 


Temperature  dependent  photorefractive  effects 

In  1988  Rytz  et  al.[9]  reported  results  of  two-wave  mixing 
experiments  performed  on  BaTiO^  crystals  in  the  temperature  range 
between  20'>  and  120®C.  They  observed  a  large  decrease  in  the 
density  of  empty  traps  as  well  as  a  significant  decrease  in  the 
response  time  at  elevated  temperatures.  These  characteristics 
are  not  in  agreement  with  the  usual  band  transport  model  of 
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Figure  17.  Calculated  power  law  dependence  of  the 
photorefractive  rise  time  versus  laser  intensity  as  a  function  of 
grating  period.  The  exponent  x  was  calculated  according  to  r  = 


photorefraction.  In  Figure  18  we  show  some  of  the  experimental 
results  obtained  by  Rytz  et  ai.[9].  In  Figure  19  we  plot  the 
behavior  expected  from  the  standard  ''Kukhtarev"  model.  The 
physical  mechanism  for  the  observed  differences  was  not  well 
understood. 

These  temperature  dependent  photoref ract ive  effects  can  be 
explained  by  secondary  centers.  At  the  intensities  of  2  to3 
W/cm^  used  by  Rytz  et  al.[9]  these  centers  are  easily  populated. 
To  demonstrate  the  effect  of  secondary  centers  the  temperature 
dependence  of  the  two-wave  mixing  gain  was  calculated  based  on 
the  secondary  center  model.  For  these  calculations  a  model  of 
one  species  of  secondary  centers  with  a  thermal  activation  energy 
of  0.5  eV  was  used.  Values  for  the  dielectric  constant  at 
different  temperatures  were  taken  from  reference  10.  The 
electro-optic  coefficient  was  assumed  to  be  proportional  to  the 
dielectric  constant  and  temperature  variations  in  the  refractive 
index  were  neglected. 

The  calculated  variations  in  the  two-wave  mixing  gain 
coefficient  with  temperature  are  shown  in  Figure  20.  These 
predictions  are  very  similar  to  the  experimental  results  of  Rytz 
et  al.[9],  and  are  consistent  with  a  decrease  in  the  effective 
empty  trap  concentration  with  an  increase  in  temperature.  This 
is  also  consistent  with  the  observed  decrease  in  magnitude  of  the 
photo-induced  absorption  with  increasing  absorption. 

Rytz  et  al.[9]  also  found  that  the  photoref ract ive  response 
time  decreased  significantly  with  temperature  increase.  the 
magnitude  of  this  effect  is  not  predicted  by  the  secondary  center 
model  if  it  is  assumed  that  the  hole  mobility  is  constant. 
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Figure  19.  Calculated  two-wave  mixing  gain  vs  grating  period  for 
different  temperatures  for  standard  single  species  "Kuhktrarev" 


model. 
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Agreement  with  experimental  results  requires  a  hole  mobility 
which  increases  with  temperature  with  an  activation  energy 
similar  to  the  thermal  ionization  energy  of  the  secondary 
centers,  characteristic  of  a  trap  dominated  mobility.  This 
suggests  that  the  secondary  centers  exhibit  a  strong  influence  on 
the  hole  mobility  and  therefore  photorefractive  response  time  in 
BaTi03 . 

VI .  SUMMARY 

In  this  effort  the  effects  of  secondary  centers  on  the 
photorefractive  properties  of  BaTi03  were  investigated  with  the 
goal  of  obtaining  enhancements  in  the  photorefractive  effects 
through  utilization  of  the  secondary  centers  by  external  pumping. 
Secondary  center  parameters  were  measured  through  photo-induced 
absorption  and  two-beam  coupling  experiments.  These  experiments 
demonstrated  that  two  species  of  secondary  centers  (denoted  by  N2 
and  N3)  were  present,  with  concentrations  of  approximately 
1.5x10^^  cm“^  and  2x10^®  cm”^.  It  was  discovered  that  at  near  IR 
wavelengths  («  800  nm)  these  secondary  centers  were  already  being 
populated  and  thus  contributing  to  and  enhancing  the 
photorefractive  effect.  Consequently,  enhancement  of  the 
photorefractive  effect  through  external  pumping  was  not  realized. 

A  numerical  method,  based  on  a  finite  difference  approach,  was 
developed  to  model  the  secondary  center  effects.  Although  the 
computation  time  required  to  implement  the  calculations  was 
prohibitive,  this  numerical  capability  has  been  successfully 
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utilized  to  calculate  photorefractive  grating  formations  in  other 
crystals  such  as  Bi3^2®^®20  which  the  calculations  are  less 
computationly  intensive. 
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